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Abstract 

We investigate the intrinsic noise of a micro-mechanical displacement detector 
based on the radio- frequency single-electron transistor (rf-SET). Using the 
noise analysis of a SET by Korotkov as our starting point, we determine 
the spectral density of the displacement noise due to the tunneling current 
shot noise. The resulting mechanical displacement noise decreases in inverse 
proportion to the increasing gate voltage. In contrast, the displacement noise 
due to the fluctuating SET island charge increases approximately linearly 
with increasing gate voltage. Taking into account both of these noise sources 
results in an optimum gate voltage value for the lowest displacement noise 
and hence best sensitivity. We show that a displacement sensitivity of about 
10~^ A and a force sensitivity of about 10~^^ N are predicted for a micron- 
sized cantilever with a realizable resonant frequency 100 MHz and quality 
factor Q ~ 10^. Such sensitivities would allow the detection of quantum 
squeezing in the mechanical motion of the micro-mechanical cantilever and 
the detection of single-spin magnetic resonance in magnetic resonance force 
microscopy (MRFM). 

PACS numbers: 85.35. Gv, 85.85. +j, 73.50.Td 
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I. INTRODUCTION 



Fast and ultra-sensitive displacement and force detection is of great interest for a broad 
range of fundamental applications. It is an essential requirement for the proposed detec- 
tion of quantum superposition statesS, quantum squeezed states!, and single-spin magnetic 
resonance in MRFmI which has also been proposed for the read-out stage of a solid-state 
NMR quantum computer.^ The common requirement of these examples is the need to per- 
form the measurement using a micro-to-nanoscale-mechanical detector at, or even below, the 
"standard quantum limit" (SQL) Axsql = {h/2muJoY^'^, in the radio-frequency range. 
Considering, for instance, a micron-sized cantilever with a resonant frequency of 100 MHz 
and quality factor Q ~ 10"^, the detection of quantum squeezed states requires a displace- 
ment sensitivity of about 10~^ A.i As another example, in detecting single-spin magnetic 
resonance, the interaction force between a nano-magnet with a targeted electron spin is 
typically on the order of about 10"^*^ N(100 aN).l If a micron-sized cantilever on which the 
nano-magnet is mounted is directly coupled at the spin precession frequency (100 MHz or 
higher), the magnetic force will cause the cantilever to oscillate at an amplitude of about 
10~^ A for a single spin. 

Following the development of the rf-SET,i^ it has been proposed to use an rf-SET- 
based displacement detector for measuring sub-angstrom motion of micron-scale and smaller 
mechanical oscillators.! In the rf-SET, 1// noise which limits the applications of the con- 
ventional dc-SET is considerably improved. A shot noise limit of 10~^ A/V Hz may be 
achievable by using an rf-SET displacement detector,! provided the coupling (gate voltage) 
between the SET and micro-mechanical oscillator is strong enough. However, the back- 
action of the SET onto the mechanical oscillator must also be taken into account:! the 
fluctuating force acting on the oscillator due to the fluctuating island charge produces an- 
other contribution to the displacement noise. The total displacement noise will be the sum 
of contributions from two sources - the shot noise in the tunneling current and the force 
noise due to the fluctuating island charge. 
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In this paper, we give a calculation of the displacement noise of the rf-SET displacement 
detector. By using the intrinsic noise analysis method of Korotkov,0 the spectral density of 
the shot noise Sj in the tunneling current and the spectral density of the force noise Sp are 
determined, and the total spectral density of the displacement noise Sx is then deduced. 
We also present two simple expressions from which good estimates of the displacement noise 
due to the shot noise and the fluctuating island charge can be obtained. Results for the 
ultimate intrinsic noise limit of the detector for different parameter values are presented and 
discussed. The displacement and force sensitivities are estimated for an optimized rf-SET 
displacement detector and we find that they meet the requirement for detecting quantum 
squeezing and single-spin magnetic resonance. 

II. INTRINSIC NOISE OF THE RF-SET DISPLACEMENT DETECTOR 

The scheme of the rf-SET displacement detector is shown in Fig. 1. The basic principle of 
the device involves locating one of the gate capacitor plates of the SET on the cantilever so 
that, for fixed gate voltage bias, a mechanical displacement is converted into a polarization 
charge fluctuation. The stray capacitance Cs of the leads contacting the SET and an inductor 
L form a tank circuit with resonant frequency ujt = (LCs)'^^"^, loaded by the SET. A 
monochromatic carrier wave is sent down the cable. At the resonant frequency the circuit 
impedance is small and the reflected power provides a measure of the SET's differential 
resistance -Rrf.0 When the gate capacitor is biased, mechanical motion of the cantilever is 
converted into differential resistance changes, hence modulating the reflected signal power. 

To simplify the algebra and make the final results concise, we shall in fact consider a 
symmetric dc-SET in the calculations with junction resistances Ri = R2 = Rj and junction 
capacitances Ci = C2 = Cj. The differences between the rf- and dc-SET noise formulae are 
inessential, with appropriate time-averages required in the formeri resulting in only small 
quantitative differences in the noise values. 

Referring to Fig. 1, when the drain-source voltage V^s is small compared to the voltage 
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e/Cs, where Cs = 2Cj + Cg is the total capacitance of the SET, and also the thermal energy 
ksT <^ eVds, the peaks of the tunneling current Ids will be well separated in gate voltage 
Vg. To derive an expression for Ids in the tunneling region, we consider only the two most 
probable island electron numbers n and n + 1, resulting in the following approximation:il 

Ids = e[bi{n) - ti{n)]a{n) + e[bi{n + 1) - ti{n + l)]a{n + 1). (1) 

The probabilities cr{n) and a{n + 1) are also given there and the tunneling rates bi{ti) from 
the bottom (top) across the ith junction of the rf-SET take the usual form [see expression 
( |X5| ) in the appendix]. 

In the previous calculation of the sensitivity of the rf-SET based detector, the shot noise 
formula Sj = 2e/ds was used. 10 This expression, which gives the maximum current noise, 
is only approximately valid when Vds is close to the tunneling threshold. Taking as our tool 
the intrinsic noise analysis of the SET by Korotkov, we obtain the following more accurate 
expression for the spectral density of the current noise Si{uj) [see appendix A for an outline 

of the deviation]: 

e^a(n)&(n + l) e''[f{n)-g{n + l)]a\n)g{n + l)~b\n + l)f{n) 
^^"^^ a{n) + b{n + l) a{n) + b{n + 1) [a{n) + bin + 1)]^ + cu^ ' 

where 

a(m) = bi{m) + t2{m), 
b{m) = b2{m) + ti(m), 
/(m) = 6i(m) - t2{m), 

gijn) = 62 (^) — ti{m). (3) 

The spectral density of the displacement noise due to the shot noise in the tunneling 
current is: 

S'M = Si{to)/{dhJdxf. (4) 

The dependence on the displacement x enters through the gate capacitance which is approx- 
imately Cg ~ ~^ x/d) for |a;| <^ d where d is the cantilever electrode-island electrode 
gap. 



The fluctuation of tlie island charge will produce a fluctuating force on the cantilever, 
therefore leading to another contribution to the displacement noise. The ultimate force noise 
is only determined by the stochastic character of the tunneling process. By using the same 
method of intrinsic noise analysis,^ the following expressions are obtained for the spectral 
density of the force noise Sp{uj) [see appendix B for an outline of the deviation]: 
AA'^e^a{n)b{n + 1) 



Sf(uj) 



a{n) + h{n + 1) 

{2Cj/Cyf{2n + 1 + 2An)2 + {ACj/Cg){2n + l)(2n + 1 + 2An) + (2n + l)^ 



[a{n) + h{n + l)f + 

where 



(5) 



A = C,{2C,-C,)/2dCl (6) 

and 

An=^-^-n--. (7) 

e 2e 2 ^ ' 

We model the cantilever as a simple harmonic oscillator, so that the spectral density of 
the displacement noise due to the fluctuating force is: 



SF{uj)/m^ 



where mefj is the motional mass of the cantilever, ujq is its resonant frequency, and Q is 
the quality factor. For a cantilever with geometry I x w x t (length x width x thickness) and 
made from material of density p and Young's modulus E, the resonant frequency of the 
lowest-order flexural mode isi 



The motional mass of the cantilever is given by rrics = plwt/A and the effective spring 
constant k^g at the cantilever's tip is kcs = rricsi^Q- The total displacement noise will be 
given as the sum of the above two contributions (§) and (|^): 

5xM = 5iM + 5|M. (10) 



In addition, the following expressions are given here only for simple estimates of the 
displacement noise at the resonant frequency of the cantilever [see appendix C for an outline 
of their derivation]: 
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X 



eV,- 



ds 



2.5kBTdC^ 



Cg{V,-V^./2) 



(11) 



C,(2Q-C,)(r,-Vds/2) 



(12) 



where ii' 40 is a constant. 



III. RESULTS AND DISCUSSION 

In this paper, we consider crystalline Si cantilevers with mass density p = 2.33 x 
10^ kg/m^ and Young's modulus Eiqq = 1.33 x 10^^ N/m^. To make the picture clear, 
we compare the displacement noise for different cantilevers. Table I lists the resonant fre- 
quencies /o and masses m (the effective mass of a cantilever is = m/4) of cantilevers with 
the same width and thickness but different lengths. The resonant frequencies /o = ujq/2t: 
are obtained from expression (^. We mainly focus our attention on the detector using the 
cantilever with geometry 1.25 x 0.6 x 0.15 in /zm, quality factor Q = 10^, resonant frequency 
/o = 117 MHz, and effective spring constant kes = 3.6 x 10^ aN/A. The standard quantum 
limit of the fundamental flexural mode of this cantilever is Axsql = 3.3 x 10~^ A, so that 
the minimum detectable force at the resonant frequency is AF = kesAx^Qi^/Q = 120 aN. 

Fig. 2 shows the dependence of the displacement noise Sx on the gate voltage Vg at 
the resonant frequency /o = 117 MHz of cantilever . The result assumes a symmetric SET 
with junction capacitance Cj = 100 aF, junction resistance Rj = 50 kQ, gate capacitance 
Cg = 50 aF, gate capacitor gap d = 0.1 fim, source-drain voltage Vds = 0.38 mV(0.6e/Cs), 
and temperature T = 74 mK(0.01e^/Cs), which are all routinely achievable. The solid 
line is obtained by using expressions (0), (H), and (|[) and the dashed line using ([TT|) . The 
minimum displacement noise occurs near the edges of the tunneling current Ids, and decreases 
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approximately in inverse proportion to increasing gate voltage. The displacement noise 
diverges in the center of tunneling current peaks and in the Coulomb blockade region where 
the tunneling current is suppressed. The result here is similar to Fig. 2 in Ref. 

The dependence of the displacement noise \fs^ on the gate voltage Vg at the resonant 
frequency /o = 117 MHz is illustrated in Fig. 3. The result assumes the same parameters as 
above. The solid line is obtained by using expressions (|) and (H), and the dashed line using 
(0). The maximum displacement noise increases linearly with the increasing gate voltage. 

Fig. 4 shows the total displacement noise \/Sx as a function of gate voltage. The tunnel- 
ing current is shown as the dashed lines. Considering the opposite dependences of Sx and 
Sx on gate voltage, we will expect an optimum gate voltage value for the lowest displace- 
ment noise and best sensitivity. But there are two problems involved in the determination. 
The first problem is that near the edges the tunneling current Ids may be too small to be 
measured accurately. The second is that because Sx and Sx change in different ways with 
increasing gate voltage Vg, it is rather involved finding the most optimum value of Vg at 
which the total noise is minimum; the positions of the minima in the displacement noise 
doublets will change relatively to the tunneling current peaks with increasing gate voltage. 
An alternative, simpler strategy is to investigate the displacement noise at the gate volt- 
age for which the tunneling current is fixed and measurable, say between 30% and 60% 
of its maximum value /ds,max- Fig- 5 shows the total displacement noise for our rf- 

SET displacement detector using different cantilever parameter values listed in Table I. The 
above method is utilized with the displacement noise determined where the tunneling cur- 
rent is half its maximum value. For the detector using a cantilever with resonant frequency 
/o = 117 MHz, the minimum displacement noise is 4.3 x 10^^ A/y/Hz, corresponding to 
an absolute displacement sensitivity of 5 x 10~^ A and force sensitivity 180 aN. Note that 
the minimum displacement noise increases and shifts to lower gate voltage for increasing 
cantilever length. As is evident from (p!2|), this is a direct consequence of the effective spring 
constant dependence: obviously a stiffer cantilever will be displaced less by the same applied 
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force. 

Fig. 6 shows the detector displacement noise for different junction capacitance values. 
Both Sx and 5*^ decrease with decreasing Cj. Because small displacements are detected by 
monitoring changes in tunneling current, taking a smaller junction capacitance serves an- 
other advantage for the experiment: with C^Vds/e kept fixed, a smaller junction capacitance 
corresponds to a larger source-drain voltage, so that the maximum tunneling current /ds,max 
increases according to expression (|C^) . 

Fig. 7 shows the displacement noise when the tunneling current is biased at various values 
between 30% and 60% of /ds,max- We see that performing detection at a smaller percentage 
of the maximum tunneling current becomes another effective way to further decrease the 
displacement noise. The lowest displacement noise we obtain is = 3.2 x 10"^ A/VEz, 
assuming Cj = 100 aF, Cg = 50 aF, and I^s = 0.57 nA, biased at 30% of its maximum value. 
The optimum gate voltage is 1.7 V in this case. 

In conclusion, approaching the SQL will become possible for a carefully engineered rf- 
SET displacement detector. A sensitivity at, or even beyond the SQL may be achievable by 
using the same parameters as in the above analysis but with the smaller junction capacitance 
Cj = 75 aF. Such a sensitivity would enable us to observe quantum squeezed states in micro- 
mechanical systems. The proposed displacement detector also meets the crucial requirement 
for advances in MRFM. The radio-frequency micro-mechanical cantilever used in our model 
allows direct coupling at the spin precession frequency, therefore allowing fast read-out, 
while the predicted force sensitivity would allow the measurement of single-spin resonance. 
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APPENDIX A: TUNNELING CURRENT SHOT NOISE 



The tunnehng current shot noise is given by expression (30) in Ref. |T0[ In our notation, 



it is: 



m,m' 

X {[t2{m - 1) - 61 (m - 1)] cr(m - 1) + [ti(m + 1) - h2{m + 1)] cr(m + 1)} , (Al) 



where 



2U = (eV2) J2 Mm) + him) + ti{m) + bi{m)] 



(A2) 



and Bfn',m is the real part of the inverse of the tridiagonal matrix {iuji — F), where / is the 
unit matrix and 



Tnm = Sn-l,ma{'m) + 5„+l,m&("^) - 6n,mC{m), 



(A3) 



with 



a{m) = hi{m) + t2{m), 
h{m) = ti{m) + b2{m), 
c(m) = 61 (m) + t2{m) + ti(m) + b2{m). 



The tunnehng rates are 



biin) = [AEr{n)/e^Rj / 
U{n) = \AE+{n)/e'Rj / 



1 _ g-AB-(n)/fcBT 



1-e" 



-AE+{n)/kBT 



where 



(A4) 



(A5) 



AE^in) = [-e/2 ± (en - QKis - CgV^)] e/Cs, 
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AE^{n) = [-e/2 t (en + (C, + Cg)V^ds - CgV^)] e/C^. 



(A6) 



A convenient form for the tunneling rates is 



61 (m) 
b2{m) 
ti(m) 
him) 



n — m + An + Vr 



ds 



RjCj] 1 — exp — (n — m + An + Vds) /T 
1 m — n — 1 — An + Vds 



RjC^ 1 - exp -(m -n-l-An + Vds) /T 
1 m — n — 1 — An — Vds 



RjCt. 1 — exp 
1 



— (m — n — 1 — An — Vds) /T' 
n — m + An — Vds 



RjCj: 1 - exp -(n - m + An - Vds) Z^" 



(A7) 



where 



An 

Vds 

T 



e 

CsVds 



2e 



— n — 



(A8) 



and n, n + 1 are the most probable island electron numbers, determined by the condition 



n < 



CgVg ^ CgVds 



< n + 1, 



which is equivalent to —0.5 < An < 0.5. 

Under the condition ksT <^ eVds (equivalently T ^ Vds), we have that ti(n), &2(n), 
6i(n + 1), t2(n + 1) ~ for small drain-source voltage, and furthermore b{m) ^ for all 
m < n and a{m) ^ for all ni > n + 1. Because the probabilities a{m) ^ for all 
m 7^ n, n + 1, only the terms m = n and ni = n + 1 need be considered in the summation 
over m in eqs. (|A1D and (|A2|) . From the form of the matrix ( [A^ ) and the above approximate 
conditions on the probabilities and tunneling rates, we then also have that only m' = n 
and m' = n + 1 need be considered in the sum over m'. The matrix elements -B„,n, -Bn,n+i, 
Bn+i,n, and B^+i^n+i contain only the tunneling rates a(n), 6(n), c(n), and c(n + 1) and it 
is valid to just consider 2x2 matrices for F and B: 
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B = Re 



luo — c[n] 



-h{n) 



^ —a{n) ioj — c{n + 1) y 
The probabilities a{n) and a{n+ 1) are given in Ref. | as 

a{n) = bin + l)/[a{n) + b{n + 1)], 
a{n + 1) = a{n)/[a{n) + b{n + 1)]. 



(A9) 



(AlO) 



Substituting ( |A^ ) and (|A10|) into (|A1|) and ( [A^ ) and neglecting all small terms, we get 
e^a{n)b{n + 1) e^[f{n) - g{n + 1)] a\n)g{n + 1) - b\n + 



+ 



a{n) + b{n + l) a{n) + b{n + l) [a{n) + b{n + 1)]'^ + u"^ 



(All) 



APPENDIX B: FORCE NOISE DUE TO THE FLUCTUATING ISLAND 

CHARGE 



The electrostatic interaction between the fixed island and cantilever electrodes of the 
gate capacitor is 



F = A[Cj{Vds - 2Vg) - ne + Q 



PI ' 



(Bl) 



where Qp is the background charge and 



A = Cgi2Cj - Cg)/2dCl 



(B2) 



Define 



Q = C,{Vds - 2Vg) = -{2n + 1 + 2An)eC,/C3, 



(B3) 



where the second equahty follows from ( |A8| ) . Neglecting Qp in ([BTD, the force noise is given 

as 

Spii^) = A' {AQ^e'Snniuj) + e^S^2^2{u) - 2Qe' [S^2^{uj) + S^A^)]} . (B4) 

Similarly to expression (27) in Ref. [1^, we have 

11 



= 4 ^ m"'Bm',mm'^(T{m). (B5) 

m,m' 

Using expressions ( |A9D , ( |A10| ), and (p5|) to obtain S'„„(c<j), S'„2„2(ct;), and [5'„2„(a;) + 5'„„2(c<j)], 
the final result for the force noise is 

4A^e^a{n)b{n + 1) 
a{n) + b{n + l) 

i2Cj/Cgf{2n + 1 + 2An)^ + (4Cj/Cg)(2n + l){2n + 1 + 2An) + {2n + 1)^ 
^ [a(n) + 6(n + 1)]2 + u;2 • ^ ^) 

Note that both Sj{uj) and 5'i?(co') have a Lorenzian spectrum. 



APPENDIX C: SIMPLE ESTIMATE FOR THE DISPLACEMENT NOISE 

A simple, approximate formula for the tunneling current i 

hs = /ds.max/ COsh^ [eCg(V"g,res " Vg) /2.5 /CfiTCs] , (CI) 

where V^,res is the value of the gate voltage at which the tunneling current is maximal: 

^.,res = ^ + (r.+ ^)^. (C2) 

Subsituting (|Cl| ) and the shot noise limit formula Sj = 2e/ds into @), we obtain 



2el 



ds,max 



2.5kBTdC^ 
CgiVg - VJ2) 



cosh^ [eCgiVg,,,, - Vg)/2.5kBTC^] 
27777777 7777777~777777 • y^'^) 



sinh^ [eCg{Vg,,,s - Vg)/2.5kBTCj,] 
Using expression (|l|) for the tunneling current with (j{n) = a{n + 1) = 0.5 and expressions 
( |A7| ) for the tunneling rates with An = gives for the maximum tunneling current: 

-^ds,max = Kis/4-Rj. (C4) 

To make a simpler estimate of the displacement noise, we might replace the dimensionless, 
parameter- independent function cosh^[a;]/ sinh^[x] by a constant say its minimum value 
which is 4. However, this function poorly approximates the correct doublet shape and it is 
found numerically that the ^ 40 provides a good fit. Thus we obtain the final result for 
the displacement noise due to the shot noise in the tunneling current: 
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ds 



2.5kBTdCs 



(C5) 



A simple approximation to the maximum force noise is obtained by setting An = in 
61) and neglecting small terms: 



^^(0) = iA'e'>R,/V^,) [Cj:{2n + l)/CgY ■ 



(C6) 



Substituting in the definition (|B2| ) for A and using 2n + 1 = Cg{2Vg — Vds)/e at maximum 
tunneling current, we obtain: 



S 



X 



Cg{2C,-Cg){Vg-V^j2 



n 2 



dCl 



(C7) 



at the resonant frequency uoq of the cantilever, which is assumed to be much smaller than 



the tunneling rates in magnitude. 
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TABLES 

TABLE I. The resonant frequency /o and mass m of crystalline Si cantilevers with different 
lengths (/) and fixed width {w = 0.6 /xm) and thickness (t = 0.15 /xm). 

I in fim 1.25 1.5 2 4 

/o in MHz 117.2 81.36 45.77 11.44 

m in pg 0.262 0.315 0.419 0.839 
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FIGURES 



FIG. 1. Scheme of the rf-SET displacement detector. 



FIG. 2. Displacement noise due to the tunneling current shot noise as a function of gate 
voltage (upper, doublet curves). The estimate of the lower limit on the displacement noise vs. 
gate voltage is also shown (lower, dashed line). The noise analysis assumes a symmetric SET at 
T = 74 mK(0.01e^/Cs) with junction capacitance Cj = 100 aF, junction resistance Rj = 50 kJ7, 
and static gate capacitance Cg = 50 aF (the gate capacitor gap is d = 0.1 /um). The drain-source 
voltage is Vds = 0.38 mV(0.6e/Cs). The displacement noise is determined at 117 MHz-the resonant 
frequency of the cantilever. 



FIG. 3. Displacement noise due to the force noise from the fluctuating island charge as a 
function of gate voltage (lower curves). The estimate of the upper limit on the displacement noise 
is also shown (upper, dashed line). The noise analysis assumes the same parameters as in Fig. 2. 
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FIG. 4. The total displacement noise as a function of gate voltage (upper, doublet curves). 
For illustrative purposes, the sample voltage range is chosen such that the force and current dis- 
placement noise contributions are comparable. The tunneling current is given by the lower, dashed 
line. 



FIG. 5. The total displacement noise as a function of gate voltage for the different cantilever 
lengths listed in Table I. The same parameter values as in Fig. 2 are assumed. The displacement 
noise at the resonant frequency of the cantilever is determined at the gate voltage for which the 
tunneling current is fixed at half its maximum value. 



FIG. 6. Total displacement noise for different junction capacitances. Except for the junction 
capacitances, the same parameter values as in Fig. 2 are assumed and the displacement noise is 
determined at the tunneling current half-maximum. 



FIG. 7. The displacement noise for different percentages of the maximum tunneling current. 
The same parameters as in Fig. 2 are used. 
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Fig. 1 Zhang 
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Fig. 3 Zhang 
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Fig. 4 Zhang 
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Fig. 5 Zhang 



22 




g' ■ ■ ■ ■ ' 

1 2 3 4 5 

Gate Voltage (V) 
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Fig. 7 Zhang 
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